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I.  OBJECTIVES 

Intermolecular  forces  govern  many  physical,  chemical,  and  biological  processes.  It  is  now 
possible  to  model  such  forces  on  computers  using  theoretical  quantum  mechanical  methods 
developed  by  our  group.  Recent  advances  in  computer  hardware  and  software  allow  inves¬ 
tigations  of  typical  molecules  found  in  explosives,  propellants,  and  solvents.  The  AASERT 
grant  supported  graduate  students  involved  in  this  research.  Their  work  helped  to  obtain 
state  of  the  art  potential  energy  hypersurfaces  for  various  molecules  relevant  for  the  defense 
mission. 


II.  RESEARCH  AND  EDUCATION  ACTIVITIES 

The  grant  period  was  7/1/97  to  6/30/00  and  the  award  amount  was  $77,689.00.  The 
grant  supported  mainly  a  graduate  student  Eric  Mas.  In  the  last  year  it  also  supported 
another  graduate  student  Garol  Murdachaev.  Both  students  were  involved  in  various  edu¬ 
cational  activities.  These  included  courses  related  to  their  research  work  taken  at  the  Uni¬ 
versity  of  Delaware,  participation  in  group  meetings  which  include  a  significant  amount  of 
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pedagogical  seminars,  and  participation  in  various  conferences  (mainly  Gordon  conferences 
on  intermolecular  forces).  Research  activities  involved  leading  edge  work  on  development 
of  intermolecular  potentials.  The  systems  studied  by  the  students,  described  in  more  detail 
in  the  next  section,  were  model  systems  enabling  development  of  methods  applied  later  by 
other  members  of  our  group  to  molecules  of  direct  interest  to  DoD,  such  as  those  involved 
in  our  SERDP  joint  project  with  ARL  [1,2].  Eric  Mas  graduated  with  a  Ph.D.  in  Physics  in 
May  1999.  He  is  now  a  postdoctoral  associate  in  the  Los  Alamos  National  Laboratory  and 
works  on  modeling  of  materials  relevant  for  the  US  defense  mission.  Clearly,  the  training 
goal  of  the  AASERT  program  could  not  have  resulted  in  a  better  outcome.  Garold  Mur- 
dachaev  continues  work  on  his  Ph.D.  research  and  is  expected  to  graduate  in  about  three 
years. 


III.  RESEARCH  ACHIEVEMENTS 

Eric  Mas  work  supported  by  AASERT  resulted  in  three  published  papers  [3-5].  The 
subject  of  these  papers  was  the  development  and  testing  of  intermolecular  potentials.  Since 
systems  of  interest  to  DoD  are  very  large  from  the  point  of  view  of  electronic  structure 
calculations,  methods  applied  to  such  systems  should  first  be  tested  on  smaller  ones.  Water 
was  chosen  as  an  excellent  test  system  since  it  is  certainly  the  most  investigated  one  and 
therefore  comparisons  with  vast  literature  data  are  possible.  Furthermore,  supercitical  water 
is  a  subject  of  significant  interest  to  DoD  and  first-principles  methods  should  play  a  partic¬ 
ularly  important  role  for  condensed  phase  under  extreme  pressures  and  temperatures.  Our 
electronic  structure  calculations  used  the  symmetry-adapted  perturbation  theory  (SAPT) 
developed  by  our  group.  Iinteraction  energies  for  about  2500  water  dimer  configurations 
have  been  calculated  using  SAPT.  Effective,  interaction-optimized  bases  were  used  leading 
to  0.3  kcal/mol  accuracy  near  the  minimum  of  the  dimer  potential.  The  computed  points 
were  then  fitted  to  two  types  of  analytic  potential  energy  surfaces:  a  site-site  form  and 
an  expansion  in  functions  dependent  on  the  vector  connecting  the  centers  of  mass  and  on 
the  Euler  angles  defining  the  orientation  of  each  monomer.  The  second  virial  coefficient 
was  calculated  from  these  surfaces  including  the  quantum  correction  and  isotopic  depen¬ 
dence,  as  well  as  the  molar  heat  capacity  at  constant  pressure.  Comparison  of  these  data 
to  experiment  shows  that  both  of  our  surfaces  are  superior  to  any  previously  available.  In 
addition,  the  three-body  nonadditivity  has  been  computed  on  a  grid  of  points  sufficient  for 
a  reduced  dimensionality  calculation  of  the  water  trimer  spectra.  The  dimer  spectrum  was 
obtained  from  the  pair  potential  through  rigorous  six-dimensional  quantum  calculations  of 
the  vibration- rotation- tunneling  (VRT)  levels.  For  the  water  trimer  three-dimensional  quan¬ 
tum  calculations  were  performed  to  reproduce  the  hydrogen  torsional  spectrum.  The  SAPT 
ab  initio  potentials  were  found  to  be  the  first  potentials  that  quantitatively  reproduce  the 
experimental  far-infrared  spectra  of  the  water  dimer.  The  three-body  interactions — together 
with  the  pair  potential — produce  an  accurate  representation  of  the  hydrogen  bond  torsional 
levels  in  the  spectrum  of  the  water  trimer.  This  work  indicates  that  first-principles  poten¬ 
tials  for  water  achieved  the  level  of  sophistication  allowing  highly  accurate  predictions  even 
for  observables  as  sensitive  as  intermolecular  spectra.  Thus,  one  can  develop  ab  initio  SAPT 
potentials  which  can  be  used  to  predict  properties  of  materials  relevant  to  the  DoD  mission 
such  as  densities  of  energetic  compounds  and  properties  of  supercitical  fluids. 
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Garold  Murdachaev  continued  the  work  on  water  with  the  aim  of  developing  a  flexible 
water  potential,  i.e.,  a  potential  depending  on  the  intramonomer  coordinates.  This  task  is 
beyond  reach  of  the  standard  ab  initio  methods  since  the  potential  energy  surface  involving 
all  degrees  of  freedom  is  12-dimensional  already  for  the  water  dimer.  However,  the  SAPT 
approach  allows  us  to  attack  this  problem  by  using  asymptotic  properties  of  the  interaction 
energy. 
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